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INTRODUCTION 


Stress  wave  attenuation  denotes  the  dissipation  of  the 
energy  of  a  propagating  wave  in  a  medium.  Every  material  is  dissi*- 
pative  in  this  respect  and  therefore,  stress  wave  attenuation  is 
a  constitutive  property  of  the  material.  The  interpretation  of 
stress  wave  attenuation  as  a  nondestructive  evaluation  (NDE)  tech¬ 
nique  in  the  sonic  range  is  as  old  as  the  skill  of  striking  a  piece 
of  glassware  or  a  metal  part  and  listening  for  the  “characteristic 
ring  of  a  good  specimen“.  The  measurement  of  stress  wave  attenua¬ 
tion  via  ultrasonic  NDE  is  a  post-World  War  II  development. 

The  primary  focus  of  ultrasonic  attenuation  NDE  techniques 
has  been  the  detection  of  overt  flaws.  The  concept  that  such 
techniques  could  be  used  to  monitor  general  structural  and  material 
degradation  and  physical  characteristics  has  received  much  less 
attention.  In  the  1950's  Truell  and  Hikata  [1]  found  that  attenua¬ 
tion  in  aluminum  increased  with  the  number  of  fatigue  cycles. 

More  recently,  Vary  et  al .  [2-7]  have  shown  that  ultrasonic  atten¬ 
uation  can  be  correlated  with  the  mechanical  and  strength  properties 
of  metals,  ceramics  and  fiber  composites.  In  testing  graphite  fiber 
polyimide  composites,  Hayford  et  al.  [8]  observed  good  correlation 
between  the  Initial  attenuation  and  the  shear  strength  as  determined 
by  the  short  beam  shear  test.  Williams  and  Doll  [9]  found  that  the 
initial  attenuation  can  be  an  effective  indicator  of  fatigue  life 
of  unidirectional  graphite  fiber  composites  subjected  to  transfiber 
compress  ion- compress  Ion  1 oad Inn, 
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Ultrasonic  attenuation  can  be  measured  experimentally  by 
either  the  pulse-echo  technique  or  the  through-transmission 
technique.  In  the  pulse-echo  technique,  one  or  two  transducers  on 
the  same  face  of  the  specimen  are  used  to  transmit  the  pulse  and 
receive  its  echo  from  a  flaw  or  the  opposite  face  [8,10,11]. 

Usually  a  short  duration  pulse  having  a  broad-band  frequency  spec¬ 
trum  Is  used.  Thus,  a  spectral  analysis  is  required  to  extract 
the  frequency  dependence  of  the  attenuation  [12,13].  The  through- 
transmission  technique  always  requires  the  use  of  separate  trans¬ 
mitting  and  receiving  transducers  placed  on  opposite  surfaces  of 
the  structure  [14-17].  Narrow-band  pulses  can  be  used  to  circum¬ 
vent  the  spectral  analysis  requirement  [17].  Through-transmission 
has  the  generally  recognized  disadvantage  that  when  thin  sections 
are  tested,  the  superposition  of  the  transmitted  signal  with  echoes, 
and  even  echoes  upon  echoes,  makes  any  interpretation  of  the  output 
very  difficult. 

The  through-transmission  technique  is  illustrated  in  Fig.  1 
where  the  transmitting  transducer  is  driven  by  a  narrow-band  fre¬ 
quency  wave-packet  electrical  pulse  which  has  a  center  frequency  co. 
Such  a  signal  is  shown  in  the  top  trace  of  Fig.  2.  A  stress  wave 
corresponding  to  the  excitation  is  introduced  into  the  specimen. 

The  stress  wave  propagates  towards  the  other  transducer  and  is 
received  as  the  output  signal.  A  typical  output  signal  is  shown 
in  the  bottom  trace  of  Fig.  2.  (The  specimen  was  a  unidirectional 
graphite  fiber  epoxy  composite  (Hercules  AS/3501-6),  1.917  cm  long 
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with  a  cross  section  of  1 .080  cm  x  I.080  cm.  The  frequency  was 
1.2  MHz.)  A  longitudinal  wave  was  propagated  in  the  X2  direction 
which  is  perpendicular  to  the  fiber  and  the  laminate  directions  as 
shown  in  Fig.  3.  A  conventional  method  for  determining  the  stress 
wave  attenuation  via  the  ultrasonic  through-transmission  technique 
involves  comparing  the  output  signals  from  two  different  lengths 
of  otherwise  identical  specimens.  The  attenuation  a(to),  which  is 
generally  frequency-dependent  can  be  calculated  as 


where  the  A.  and  the  £.  are  the  output  signal  amplitudes  and  the 
lengths  of  the  respective  specimens,  and  Zn  is  the  natural  logarithm. 

If  the  time  required  for  the  input  pulse  to  traverse  the 
length  of  the  specimen  is  much  less  than  the  pulse  width,  the  out¬ 
put  signal  will  contain  overlapping  echo  reflections.  Such  an 
example  is  shown  in  the  bottom  trace  in  Fig.  A.  (The  specimen  was 
the  same  as  that  tested  in  Fig.  2  except  that  the  length  was  0.401  cm 
and  the  frequency  was  0.5  MHz.)  The  initial  portion  of  the  output 
signal  is  complicated  by  both  the  transient  response  of  the  receiver 
transducer  and  the  superposition  of  echo  reflections  while  the  latter 
portion  of  the  output  signal  contains  several  echo  reflections. 

Thus,  the  output  signal  is  difficult  to  interpret  in  terms  of  the 
material  attenuation.  It  is,  however,  very  interesting  to  note  that 
the  output  signal  settles  to  a  steady-state  oscillation.  This 
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report  describes  the  use  of  this  easily  measured  steady-state 
output  amplitude  to  determine  the  attenuation  of  the  specimen. 


STEADY-STATE  OUTPUT  AMPLITUDE  ANALYSIS 


In  the  through-transmission  technique  illustrated  in 
Fig.  1,  a  sinusoidal  electrical  signal  of  voltage  amplitude  Vj 
is  sent  into  the  transmitting  transducer  which  introduces  a  propa¬ 
gating  stress  wave  aj(x,t)  given  by 

a^(x,t)  =  V jF^ (w)e' (2) 

into  a  nondissipative  specimen.  The  time-varying  character  of  the 
signal  is  represented  by  the  complex  notation  where  i  =  and  only 
the  real  part  of  eqn.  (3)  should  be  considered.  The  coordinate  x 
is  positive  when  pointing  away  from  the  transmitting  transducer  and 
towards  the  receiving  transducer  (see  Fig.  1),  and  t  is  time.  The 
proportionality  constant  F^ (u)  is  a  function  of  co  and  is  formally 
the  transduction  ratio  for  the  electromechanical  gyrating  system 
[l8]  corresponding  to  the  transformation  of  an  input  electrical 
voltage  to  a  mechanical  stress  wave  in  the  specimen.  F^(a)) 
includes  both  the  characteristics  of  the  transmitting  transducer 
and  the  stress  wave  transmissibi 1 i ty  of  the  transducer-specimen 
interface.  The  parameter  F^  (o))  has  the  dimensions  of  stress  per 
unit  voltage.  The  wave  number  k  is  defined  by 


where  the  wavelength  X  is  given  by 

C^(a)) 

X 


(3) 


(4) 
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and  where  Cp(w)  is  the  phase  velocity  and  f  )  is  the  fre¬ 

quency  in  Hertz.  For  wave  propagation  in  a  dissipative  medium,  the 
stress  wave  attenuation  may  be  expressed  in  terms  of  the  attenuation 
parameter  a(a)),  and  so  eqn.  (2)  should  be  replaced  by 

(x,t)  =  v.F^e  e  '  .  (5) 


The  primary  stress  wave  is  defined  as  the  initially  propagated 
stress  wave  between  x  =  0  and  x  =  Z.  The  primary  stress  wave  at 
the  end  of  the  specimen  x  =  Z  is  given  by 

The  electrical  output  signal  from  the  receiving  transducer  corre¬ 
sponding  to  this  primary  wave  is 


Vj  (t) 


_  _  -a£  i  (k£  -  oot  -  0) 

v^FiF^e  e 


(7) 


where  0(aj)  Is  a  phase  shift  (O  £  9  <  27r) ,  and  the  proportionality 
constant  F^((jo)  is  the  transduction  ratio  corresponding  to  the  trans¬ 
formation  of  a  stress  wave  in  the  specimen  to  an  output  electrical 
voltage  and  Includes  both  the  characteristics  of  the  receiving 
transducer  and  the  stress  wave  transmi ss i b 1 1  I ty  of  the  transducer- 
specimen  interface.  The  parameter  has  the  dimensions  of 

voltage  per  unit  stress,  and  thus  the  product  Fj(a3)F2(co)  is 
dimensionless. 

A  portion  of  the  stress  wave  in  eqn.  (6)  reflects  from  the 
end  of  the  specimen  with  a  reflection  coef  f  ic  i  ent  R  (co) ,  (0  <  R  <  1), 
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and  is  given  by 


a2(x,t) 


-a{2l  -  x)  i[k(2£  -  x)  -  u)t  -  (J)] 
V . r , Ke  e 

I  1 


(8) 


where  ({)(u),(0  <  (p  <  Zir) ,  is  a  phase  shift  due  to  the  boundary 


reflection. 

The  analysis  of  this  process  can  be  continued  indefinitely 
by  using  the  following  rules:  (1)  a  stress  wave  reflected  from  an 
end  of  the  specimen  can  be  obtained  from  the  incident  stress  wave 
by  Imposing  a  phase  shift  of  <j>  and  by  multiplying  the  wave  ampli¬ 
tude  by  the  reflection  coefficient  R;  and  (2)  an  electrical  output 
signal  can  be  obtained  from  the  stress  wave  arriving  at  the  receiver 
end  of  the  specimen  by  imposing  a  phase  shift  of  9  and  by  multi¬ 
plying  the  wave  amplitude  by  the  parameter  F^.  This  process  is 
illustrated  in  Fig.  5* 

The  steady-state  portion  of  the  output  signal  of  the 

receiving  transducer  v*  consists  of  the  transduced  portion  of 
^  0  ,ss 

the  primary  stress  wave  and  the  transduced  portions  of  all  its 
reflections,  and  is  given  by 


^  ^  ^  -a£  i(k£  -  cot  -  0)  ^  P  P  2  -3a£^i(3l<£ 

V*  =  v.F,F.,e  e  '  +  v.F,F,R  e  e 

o,ss  I  1  2  I  I  2 


V  „2j  -(2j  +  Doe  i[(2j  +  l)k£  -  wt  -  2j(() 

''i'■^2  Z 
j  =  0 


-  cot  -  2^)  -  9) 


-  0]  .  (9) 
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If  the  signal  represented  by  eqn.  (9)  is  phase-shifted  by  an  amount 
(k£  -  0),  the  phase-shifted  steady-state  output  becomes 

V  =  v.F.F,  y  *  l)rf  i[2J(kf  -  «  -  ml  ,  ^ 

o,ss  I  1  2  Z_( 

j  =  0 

No  generality  is  lost  by  using  eqn.  (10)  instead  of  eqn.  (9)  because 
the  exact  phase  information  of  the  output  signal  is  not  used.  By 
def in ing 

(p‘(o})  =  kl  -  (P  ,  (11) 
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The  summations  in  eqn.  (14)  can  be  summed  exactly  to  [19]" 


V  =  v.F.F.e 
o,ss  I  1  2 


(Xl[ 

’  "I 

1  cosoot 

1 

L 

sinh(2a£  -  2£nR) 
cosh(2a£  -  2lnR)  -  cos(2(})') 


-  1 


.  ^  ,,  .  _ sin(2(t) ') _ _ 

+  s I nwt  2  cosh(2o6^  -  2£nR)  -  cos (24)') 


which  can  be  rewritten  as 


''i'^lV 


-CLl 


0  ,ss  2{cosh[2(a^  -  -tnR)]  -  cos(2(j)')> 

-  coswtcos (2(|) ' )  +  s intots in (241 ' ) 


Z{al  -  £nR) 
e  coswt 


(15) 


Extrema  of  v  can  be  obtained  by  using 
0  9  ss 


dv 


0  ,ss  _ 
dt 


0  . 


(16) 


The  useful 


equations  are 


2 


e  ’^\in{n6)*j 


_ sin6 

coshy  -  cos6 


n  =  0 


Y  >  0 


1  +  2 


2 
=  0 


-ny  .  sinhy 

®  ^0=("-5)=^03hY  -■  cos^ 


y  >  0 


Note  that  because  (2at  -  2£.nR)  >  0,  these  relations  are  appropriate 
for  the  desired  summations. 
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Substitution  of  eqn,  (15)  into  eqn.  (16)  gives 


s  incot 
coscot 


tancot  = 


sin(2(})' ) 
2(al  -  £nR) 


-  cos(2(f>') 


(17) 


Noting  that 


cos^6 


-  1  +  tan  0,  eqn.  (17)  can  be  rewritten  as 


COS(Ot  = 


1  + 


s  i  n  (2(fi ' ) 
2(aZ  -  ^nR) 


(18) 


-  cos(2(J)') 


Substitution  of  eqns.  (17)  and  (l8)  into  eqn.  (15)  and  manipulation 
of  the  resulting  expression  give  the  ampl itude  of  the  steady-state 
output  as 

V  F  F  e'°^  -  £nR)  _  2cos(2(})')e^^°^  “  +  1 

max  V  = - - - - - — - -  .  (]9) 

2{cosh[2(a£  -  £nR)]  -  cos(2(;)')} 

A  strategy  for  using  eqn.  (19)  to  determine  the  attenuation  of 
materials  from  the  measured  through-transm i ss ion  steady-state  output 
amplitude  will  be  illustrated  in  the  next  section. 
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EXAMPLE  OF  APPLICATION 


If  the  parameters  (ti))F2((i)) ,  R(tj)  and  (}>(('^)  s^e  known, 
the  attenuation  a(u))  can  be  determined  from  eqn.  (19)  and  the 
experimentally  measured  steady-state  output  amplitude.  Note  also 
that  the  phase  velocity  Cp(w)  is  required  for  calculating  the  wave 
number  k. 

In  general,  the  required  parameters  can  be  measured  a  priori 
by  testing  different  lengths  of  otherwise  identical  specimens.  For 
each  frequency  of  interest  the  phase  velocity,  Cp(to)  ,  and  the 
parameters  a{u)),  F^  (o)) F2((d)  ,  R(w)  and  <p(a))  can  be  found  in  that 
ord^r  based  on  the  experimental  data.  Then,  Cp{w),  F^(u3)F2(w)»  R((^) 
and  (j)(u)  can  be  assumed  constant  for  future  tests  on  the  material 
where  the  specimens  have  the  same  cross  section  but  any  length. 

An  example  demonstrating  these  calculations  for  some  graphite 
fiber  epoxy  composite  samples  follows. 

The  material  was  a  unidirectional  graphite  fiber  epoxy 
composite  (Hercules  AS/3501-6).  The  specimens  were  rods  with  a 
square  cross  section,  I.080  cm  (0.425  in)  on  a  side  as  shown  in 
Fig.  3.  The  lengths  of  the  specimens  were  0.401  cm  (O.I58  in), 

0.662  cm  (0.261  in),  I.O86  cm  (0.42?  in),  1.2/8  cm  (0.503  in), 

1.917  cm  (0.755  in),  2.545  cm  (1.002  in),  3-797  cm  (1.495  in)  and 
5.093  cm  (2.005  in).  The  direction  of  longitudinal  wave  propagation 
was  perpendicular  to  both  the  fiber  and  the  laminate  directions 
and  is  indicated  in  Fig.  3  as  the  X2  direction.  The  experimental 
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system  is  shown  in  Fig.  6.  The  transducers  were  Acoustic  Emission 
Technology  (AET)  FC-500  transducers  and  the  couplant  was  AET  SC-6. 
The  peak-to-peak  input  voltage  to  the  transmitting  transducer  was 
280  volts,  A  pressure  of  0.3  MN/m  (44  psi)  was  applied  to  the 
transducer-specimen  Interface  through  the  spring  scales.  As  shown 
in  [17],  this  pressure  exceeded  the  ‘'saturation  pressure",  which  is 
defined  as  the  minimum  interface  pressure  that  results  in  the  max¬ 
imum  output  signal  amplitude,  all  other  parameters  being  held 
constant.  The  frequency  range  was  0.3  3MHz. 

The  wave  group  velocity  C  was  determined  by 


(20) 


where  t.  are  the  wave  packet  transit  times  in  the  respective 
specimen  lengths  .  The  experimentally  measured  group  velocity 
is  independent  of  frequency  and  is  equal  to  2.82  x  10^  cm/sec 
(1.11  X  10^  in/sec).  Because  the  group  velocity  is  independent 
of  frequency,  the  composite  material  is  nond i spers i ve,  and  thus 
the  group  velocity  is  equal  to  the  phase  velocity.  Using  suffi¬ 
ciently  long  specimens  in  order  to  achieve  nonoverlapping  pulses 
such  as  depicted  in  Fig.  2,  the  attenuation  a(a3)  of  the  material 
was  calculated  from  the  amplitudes  of  the  primary  stress  wave 
output  signals  using  eqn.  (1)  and  is  plotted  in  Fig,  7. 

The  parameter  Fj(a))F2(ca)  can  be  obtained  from  the  amplitude 
of  the  primary  stress  wave  output  signal  also.  Considering  only 
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the  amplitude  of  the  signal  in  eqn.  (7)  and  rearranging  eqn.  (7) 
give 

F,(a3)F,(w)  .  (21) 

Using  the  attenuation  plot  in  Fig.  7  and  the  experimental  values 
of  v^  and  v.,  the  resulting  F^(co)F2(aj)  can  be  plotted  as  shown  in 
Fig.  8. 

In  order  to  determine  the  reflection  coefficient  R(a)),  it 
is  necessary  to  measure  the  first  reflection  amplitude 
from  the  output  signal  when  the  second  reflection  amplitude  is 
negligible.  Such  a  situation  is  illustrated  in  Fig.  9  where  the 
length  of  the  specimen  was  1.086  cm.  The  amplitude  of  the  first 


reflection 


reflect 


as  indicated  In  Fig. 


_  _  -2  -3al 
V  XI  =  v.F.F.R  e 
reflect  i  1  2 


Rearranging  eqn.  (22)  gives  the  reflection  coefficient 


Again,  note  that  in  accordance  with  energy  conservation,  the  reflec 


tion  coefficient  is  subject  to  the  restriction  0  _<  R  ^  1  .  Using 


Figs.  7  and  8  and  the  experimental  values  of  and  v.,  the 

resulting  R((jo)  can  be  computed  and  is  shown  in  Fig.  10. 
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The  phase  shift  parameter  (j)(aj)  can  be  determined  from 
eqn.  (19)  based  on  the  steady-state  output  amplitude.  Good 
accuracy  can  be  achieved  if  the  steady-state  output  consists  of 
sizable  reflections  such  that  the  primary  stress  wave  does  not 
dominate  the  output  signal.  Using  the  curves  in  Figs.  7,  8  and  10, 
a  value  of  cos(2(f)’)  can  be  obtained  from  the  steady-state  output 
amplitude  and  eqn.  (19)  by  trial  and  error.  The  phase  shift  c|)  can 
then  be  computed  from  the  arc  cosine  and  eqn.  (11).  However, 
because 

cos  ^(2(j)’)  =  cos  ^  (27T  -  2<()’)  ,  (24) 

these  calculations  result  in  two  different  values  for  ({)’  and  thus 
Substitution  of  the  two  values  of  (()  into  eqn.  (19)  to  compute  the 
steady-state  output  amplitude  for  a  different  length  specimen  and 
comparison  of  the  computed  value  with  the  experimental  value  will 
select  the  correct  (}).  The  resulting  phase  shift  parameter  (f)(a))  is 
plotted  i n  Fig .  11. 

With  the  use  of  Figs.  8,  10  and  11,  it  is  possible  to 
compute  the  attenuation  a(a))  from  eqn.  (19)  by  measuring  the 
steady-state  output  amplitude  from  any  length  (thickness)  specimen 
of. the  material  assuming  that  the  same  equipment,  specimen  cross 
section  and  direction  of  longitudinal  wave  propagation  are  maintained. 
Of  course,  the  assumption  is  that  the  phase  velocity  and  the  param¬ 
eters  F^(a3)F2((Ji)),  R(oo)  and  (t)(a))  remain  invariant  with  specimen 
substitution.  It  has  been  shown  [9]  that  the  phase  velocity  when 


compared  with  the  attenuation  is  relatively  insensitive  to  the 
fatigue  state  of  Hercules  AS/3501-6.  Further,  if  the  surface 
condition  is  consistent  from  specimen  to  specimen,  the  parameters 
(a))F2(ti)) ,  R(a))  and  <{i((jo)  are  expected  to  remain  significantly 
invariant . 
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CONCLUSIONS 


The  steady-state  output  amplitude  of  an  ultrasonic  through- 
transmission  measurement  has  been  analyzed  and  the  result  has  been 
given  in  closed  form.  Assuming  that  a  set  of  similar  specimens  of 
the  same  cross  section  can  be  tested  a  priori  in  the  same  wave 
mode  in  order  to  tabulate  values  of  the  product  of  the  input  and 
output  transduction  ratios  (F^  (a))F2(w) ) ,  the  specimen-transducer 
reflection  coefficient  (R(a))),  the  specimen-transducer  phase  shift 
parameter  (4>(w)),  and  the  material  phase  velocity  (Cp((*))),  the  atten¬ 
uation  a((jj)  for  future  specimens  can  be  determined  from  the  steady- 
state  output  amplitude  and  eqn,  (19). 

The  steady-state  output  amplitude  has  the  advantage  of 
being  easily  measurable.  Because  echo  reflections  are  actually 
accounted  for,  there  is  no  need  to  avoid  overlapping  echoes  in 
the  output  signal  and  thus  there  is  no  apparent  limitation  on  the 
minimum  length  or  thickness  of  the  specimen.  An  additional  advantage 
of  using  the  steady-state  output  amplitude  is  that  once  Figs.  8,  10 
and  11  have  been  obtained,  further  testing  does  not  require  the 
more  specialized  pulse  oscillator  or  a  spectrum  analyzer  because 
a  simple  harmonic  signal  generator  is  sufficient. 
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Fig,  1  Schematic  of  specimen  and  transducers  in 
ultrasonic  through-transmission  testing. 
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Top  trace:  Input  signal  to  the  transmitting  transducer 
Scale  =  5v/Large  Div. 

(The  signal  was  attenuated  29dB  before 
displaying.) 

Bottom  trace:  Output  signal  from  the  receiving  transducer 

Scale  =  O.lv/Large  Div. 

Time  scale:  2ps/Large  Div. 

Frequency :  1.2  MHz 


Fig.  2  Typical  ultrasonic  through-transmission  input 
and  output  signals  with  no  echo  superposition. 
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Laminate  Direction 
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Top  trace:  Input  signal  to  the  transmitting  transducer 
Scale  =  5v/Large  Div. 

(The  signal  was  attenuated  29dB  before 
displaying.) 

Bottom  trace:  Output  signal  from  the  receiving  transducer 

Scale  =  0.2v/Large  Div. 

Time  scale:  5Ms/Large  Div. 

Frequency:  0.5  MHz 


Fig.  k  Typical  signals  in  ultrasonic  through-transm i ss ion 
testing  showing  the  complicated  output  waveform  due 
to  transients  and  multiple  reflecticms  and  showing 
the  eventual  steady-state  amplitude. 
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Schematic  of  experimental  system. 
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Top  trace:  Input  signal  to  the  transmitting  transducer 
Scale  =  5v/Large  Div. 

(The  signal  was  attenuated  29dB  before 
displaying.) 

Bottom  trace:  Output  signal  from  the  receiving  transducer 

Scale  =  0.2v/Large  Div. 

Time  scale:  5us/Large  Div. 

Frequency :  1.2  MHz 


Fig.  9  Typical  ultrasonic  through-transmission  output 
signal  illustrating  first  reflection. 
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